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ASaTRACT 


A theoretical  Investigation  has  been  made  to  determine  the  variation  of  lift 
with  time  on  a slender  elliptic  cylinder  oscillating  In  pitch.  The  development 
of  the  unsteady  two-dimensional  laminar  boundary  layer  over  the  surfaces  of  the 
pitching  ellipse  Is  calculated  and  the  periodic  variation  of  lift  Is  determined 
by  matching  the  rate  at  which  boundary  layer  developed  vortlclty  Is  shed  Into  the 
wake  with  the  time  rate  of  change  of  circulation  about  the  ellipse.  The  effects 
of  mean  angle-of-attack  and  oscillation  frequencies  on  lift  hysteresis  loops  are 
determined.  It  Is  shown  that  the  hysteresis  loops  change  direction  for  mean 
angles  of  attack  greater  than  that  corresponding  to  maximum  steady  state  lift. 


2, 


TABLE  OF  CONTENTS 


List  of  Figures 1 

List  of  Tables 11 

Abstract  1 

Introduction  1 

Analysis  ^ 

Potential  Flow 4 

Governing  Equations  6 

Circulation  Criterion  14 

Solution  Procedure  16 

Discussion  of  Results  20 

Summary 26 

References 28 

Nomenclature 29 

£ 


LIST  OF  FIGURES 


Figure  1 
Figure  2 
Figure  3 
Figure  4 

Figure  5 
Figure  6 

Figure  7 


Basic  Coordinate  System  and  Body  Orientation 

Steady-State  Life  Curve  for  the  Elliptic  Cylinder 

Unsteady  Lift  Hysteresis  for  Several  Values  of  Oq  with  u ■ 0.10 

Unsteady  Lift  Hysteresis  for  a Vertal  23010-1.85  Airfoil  at 
Constant  Oscillation  Frequency  (Reference  2) 

Variation  of  with  Frequency  at  ■ <5* 

Movement  of  Upper  and  Lower  Surface  Separation  Points  for 

a =9* 

0 ^ 

Contrasting  Hysteresis  Loop  Directions  for  (u  *■  0.10. 


1 


MhiiiMfaiiimiii  irfiinfr"  rr 


< iJL 


i 


LIFT  HYSTERESIS  OF  AN  OSCILLATING 
SLENDER  ELLIPSE 


ABSTRACT 

A theoretical  investigation  has  been  made  to  determine  the  variation 
of  lift  with  time  on  a slender  elliptic  cylinder  oscillating  In  pitch.  The 
development  of  the  unsteady  two-dimensional  laminar  boundary  layer  over  the 
surfaces  of  the  pitching  ellipse  Is  calculated  and  the  periodic  variation 
of  lift  Is  determined  by  matching  the  rate  at  which  boundary  layer  developed 
vorticlty  Is  shed  Into  the  wake  with  the  time  rate  of  change  of  circulation 
about  the  ellipse.  The  effects  of  mean  angle-of- attack  and  oscillation 
frequencies  on  lift  hysteresis  loops  are  determined.  It  Is  shown  that  the 
hysteresis  loops  change  direction  for  mean  angles  of  attack  greater  than  that 
corresponding  to  maximum  steady  state  lift. 


INTRODUCTION 
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In  recent  years,  several  practical  problems  encountered  In  the  design 

of  helicopter  rotors  and  axial  compressors  have  spurred  Interest  In  the 

effect  of  tmsteady  stall  on  the  lift  of  blades  undergoing  periodic  changes 

In  angle-of-attack.  The  variation  of  lift  with  angle-of-attack  for  unsteady 

flow  Is  vastly  different  In  the  stall  region  from  the  familiar  static  lift 

curve.  In  the  steady  case  there  Is  a sharp  decrease  In  lift  at  ragles-of- 

attack  greater  than  the  stall  angle,  whereas  In  the  unsteady  oscillatory 

case  tests  by  Halfman  et  al^^^  Indicate  that  the  maximum  dynamic  lift  co- 

(2) 

efficient  Is  greater  than  the  static  C^^  . Lllva  et  al  experimentally 

max 

obtained  hysteresis  loops  In  the  vicinity  of  the  steady  stall  angle  of  attack 
on  two-dimensional  airfoils  sinusoidally  oscillating  about  a mean  angle-of 


•»'l  I 


attack.  This  study  Indicates  that  lift  hysteresis  Is  strongly  Influenced 


by  frequency,  oscillation  amplitude,  and  mean  angle-of-attack. 

Several  Investigations  have  been  performed  which  are  Invlscld  In  basic 

(3) 

formulation  but  Include  empirical  estimates  of  the  viscous  effects.  Ham 

developed  a theory  for  two-dimensional  dynamic  stall  by  modeling  the  process 

(4) 

with  discrete  vortices  shed  from  the  leading  edge.  Ericsson  and  Reding 
developed  a quasi-steady  analytical  method  In  which  airfoil  static  aerodynamic 
data  may  be  used  to  predict  the  unsteady  variations  In  lift  and  moment.  While 
methods  such  as  these  may  be  useful  In  making  predictions  of  tinsteady  aero- 
d3mamlc  effects,  they  do  not  go  to  the  root  of  the  problem.  Basically,  stall 
Is  a phenomenon  which  Is  directly  related  to  boundary  layer  separation  and 
any  complete  analysis  of  stall,  either  steady  or  unsteady,  must  Include  the 
effect  of  boundary  layer  separation. 

Several  Investigators  have  modeled  unsteady  stall  Including  the  effects 
of  boundary  layer  development.  Crlml^^^  analyzed  the  flow  about  an  airfoil 
undergoing  unsteady  motion  by  modeling  the  basic  flow  elements  and  considering 
Interactions  between  the  Invlscld  and  viscous  flows  by  Iterative  techniques. 
His  work  seems  limited  by  a linearized  description  of  the  potential  flow  which 
Is  solved  numerically  for  each  time  step  by  considering  a distribution  of 
source  and  sink  singularities.  Patay^^^  used  a boundary  layer  analysis  to 
study  the  effect  of  leading  edge  separation  on  the  delay  of  dynamic  stall  for 
pitching  Joukowskl  airfoils.  This  study  employs  a momentum  Integral  procedure 
but  solves  only  a quasi-steady  description  of  the  flow  process.  Moore^^^ 
considered  the  problem  of  predicting  the  lift  on  a slender  elliptic  cylinder 
In  a free  stream  oscillating  In  angle-of-attack.  In  addition  to  considering 
a stationary  cylinder  In  an  oscillating  stream,  Moore  simplified  his  analysis 
by  limiting  the  problem  to  the  case  of  maxlmtnn  static  lift  and  by  considering 
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only  small  rates  of  change  of  Incidence.  By  restricting  his  analysis  to 
maximum  steady  state  lift,  Moore  eliminates  the  quasi-steady  variation  In 
circulation.  In  the  lift  versus  angle-of-attack  curve,  Moore  obtains  a 
hysteresis  loop  whose  direction  Is  counterclockwise.  In  contrast  to  the 
clockwise  direction  of  hysteresis  loops  obtained  experimentally  on  cylinders 
of  alrfoli  ~ross-sectlon.  Sears'  ’ reviewed  Moore's  work  and  suggested 
that  further  Investigations  of  this  type  be  made.  Sears  discussed  the 
effects  of  Moore's  assumptions  and  stated  that  "...  the  dominant  effect  in 
this  hysteresis  may  arise  from  quasi-steady  terms  that  Moore  eliminated  by 
considering  small  perturbations  from  maximum  lift." 

The  present  study,  which  Is  not  restricted  to  maximum  steady  state  lift 
and  Includes  the  quasi-steady  variations  In  circulation,  examines  the  unsteady 
viscous  flow  about  a 6:1  elliptic  cylinder  oscillating  periodically  In  a sta- 
tionary uniform  stream.  The  major  reason  for  choosing  this  type  body  is  that 
the  Invlscld  flow  about  an  elliptic  cylinder  is  analytic  and  easily  expressed 
in  closed  form.  The  analysis  Includes  representation  of  the  unsteady  velocity 
components  as  a steady-state  term  plus  the  first  term  of  a power  series  expan- 
sion as  well  as  transformation  of  the  coordinate  system  such  that  the  separation 
point  appears  steady  In  the  transformed  coordinate.  The  differential  equations 
resulting  from  the  boundary  layer  formulation  were  solved  using  an  unsteady 
momentum  Integral  technique.  The  unsteady  circulation  as  a function  of  the 
periodically  varying  angle-of-attack  was  determined  by  iteration  using  the 
revised  Taylor-Howarth  criterion  for  relating  the  rate  at  which  boundary 
layer  developed  vorticlty  is  shed  into  the  wake  to  the  time  rate  of  change 
of  circulation. 

The  results,  which  are  presented  for  a range  of  mean  Incidence  angles 
and  oscillation  frequencies,  describe  the  behavior  of  the  unsteady  boundary 
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layer  and  Illustrate  dynamic  stall.  Including  lift  hysteresis. 

ANALYSIS 

The  coordinate  system  and  body  orientation  used  in  the  present  study 
is  presented  in  Figure  1.  The  x coordinate,  which  is  measured  from  the 
stagnation  point  for  steady  flow  at  the  mean  angle-of-attack,  is  positive 
on  the  upper  surface  and  negative  on  the  lower  surface.  The  angle  n is 
measured  positive  counterclockwise  from  the  trailing  edge  in  a body  fixed 
system.  The  angle-of-attack,  a(t),  is  assumed  positive  when  the  leading 
edge  is  pitched  up  relative  to  the  direction  of  the  stream. 

Potential  Flow 

The  potential  flow  over  an  oscillating  elliptic  cylinder  in  a free-stream 
was  obtained  from  a form  of  the  complex  velocity  potential  similar 

(9) 

to  that  given  by  Kochin  . The  complex  potential  used  in  the  present  analysis 
and  given  by 

U - t")  T' 

w * - [b  cos  a + i a sin  a] ^^7 log  ? (1) 

was  obtained  by  mapping  the  flow  about  a circular  cylinder  onto  that  around 
an  elliptic  cylinder.  With  standard  potential  flow  procedures  and  by  letting 

u'  j, , Q'L 

“6  ” “U  ’ ^ “ 2iTUL  ’ “z  “ ~ ’ 

the  Invlscld  flow  about  an  oscillating  infinite  elliptic  cylinder  was  found 
to  be: 

Ug(x,t)  » — [(1  + B)  sin  (n  - o(t))  + r(t)] 

(2) 

1/2 

where  R ■ (sin^  n + B^  cos^  n) 
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The  pocentlal  flow  for  a stationary  elliptic  cylinder  in  an  oscillatory  stream 
may  be  obtained  from  equation  (2)  by  omitting  the  term  containing  The 

resulting  expression  for  Ug(x,t)  is  then  Identical  with  that  studied  by 
Moore^^^  in  his  analysis  of  dynamic  stall.  The  angle-of-attack  and  the 
circulation  are  each  fractions  of  time  and  are  assumed  to  have  the  following 


o(t)  “ + Aa  sin  wt 

r(t)  = Aa  cos  (ot  + Ao  sin  ut 


(3) 

(4) 


I 
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where  Au  is  a constant  perturbation  of  the  Incidence  angle.  The  pitch  rate, 
is  given  by 


Q s . d(t)  * - 0)  Aa  cos  (ot 
2 

The  expression  for  r(t)  is  the  sum  of  sine  and  cosine  components,  thus 
permitting  phase  lead  and  lag  in  the  final  expression  for  circulaclon. 
Substituting  (3)  and  (4)  Into  (2)  and  discarding  terms  of  order  (Aa)^  or 
higher  yields: 

u,(x,t)  = u,.(x)  + u.(x)Ao  cos  (Ot  + u,(x)Aa  sin  (ot 
<^0  ^10  ^11 


where : 


(1  + e)sin(n  - a^)  + 

R 


(5) 


(6) 

(7) 


- (og  - ^-cos  2n 
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- (1  + 0)cos(n  - a^) 


(8) 


(9) 


"11 

The  boundary  layer  streamwlse  coordinate,  x,  is  related  to  the  angular 
measure,  n,  by 


dx/dn  ■ - R(ti)/2. 


6 


Governing  Equations 

The  non-dlmenslonal , unsteady,  laminar  boundary  layer  equations  for 
Incompressible  flow  In  a moving  coordinate  system,  provided  frequency  Is 
sufficiently  small,  are  given  by: 

+|i-  0 

3x  3y 


(10) 


3u  , . 3u  ^*^6  , ^*^6  3^u 

3t  3x  3y  " 3t  '*6  3x  3y^ 


(11) 


with  boundary  conditions 


u(x,0,t)  « 0 , v(x,0,t)  « 0 

11m  u(x,y,t)  * Ug(x,t).  (12) 

y-H» 

Physically,  the  points  of  boimdary  layer  separation  on  the  upper  and 
lower  surfaces  move  back  and  forth  with  time  as  the  body  oscillates.  The 
present  analysis  accounts  for  this  effect  by  transforming  the  boundary  layer 
equations  (10)  and  (11)  with  a scaling  of  the  x-coordlnate  so  that  separation 
always  occurs  at  a constant  value  of  the  streamwlse  coordinate  In  the  new 
system.  Before  development  of  this  transformation,  a brief  discussion  of  the 
nature  of  unsteady  separation  seems  necessary. 

Moore^^^\  Rott^^^\  and  Sears^®\  all  noted  that  vanishing  of  the 
wall  shear  does  not  necessarily  Imply  separation  In  unsteady  boundary  layer 
flow.  Instead,  they  Independently  postulated  a model  (known  as  the  HRS 
model)  In  which  the  unsteady  separation  point  Is  characterized  by  the  simul- 
taneous vanishing  of  the  shear  ai^d  the  velocity  at  some  point,  usually  away 
from  the  wall,  as  seen  by  an  observer  moving  with  the  separation  point.  The 
exact  Implications  of  this  definition  of  unsteady  separation  on  the  analytical 
prediction  of  separation  point  location  Is  not  known.  However,  Crlml^^^  Indicates 
that  for  typical  helicopter  rotor  Reynolds  Number  and  dimensionless 


i 


pitch  rate,  the  difference  between  streanwlse  location  of  the  separation 
point  as  given  by  the  MRS  model  and  the  location  of  the  classical  point 
of  vanishing  wall  shear  Is  approximately  0.6  percent  of  the  airfoil  chord. 
This  result.  In  addition  to  the  physical  reality  that  as  pitch  rate 
approaches  zero,  the  streamwlse  location  of  unsteady  and  steady  separation 
points  must  coincide,  leads  to  the  adoption  of  vanishing  wall  shear  as 
adequately  defining  the  location  of  separation.  Even  though  In  the  strict 
sense  this  definition  is  Inaccurate,  for  the  lower  range  of  pitch  rates 
under  consideration  here  this  criterion  is  well  within  the  accuracy  of  an 
analysis  to  order  Aa. 

To  scale  the  x-coordlnate,  the  new  streamwlse  coordinate, 


Is  given  by 


5 = xq(t) 


where  q(t)  Is  some  fimctlon  to  be  determined  which  adjusts  the  streamwlse 
distance  so  that  separation  always  occurs  at  a constant  value  of  C*  With 
this  scaling,  the  governing  equations  (10),  (11),  and  boundary  conditions 


(12)  become: 


3u  , 8v  - 


3u  . ^ X X ^ ^6  ^6  *^6  3^u 

— + ■“  — + uo  — + V — “ + + au  ■■  + — ^ 

3t  q 3C  ^ 35  3y  q 35  3t  ^ 6 35  ^ 


u(5»0,t)  - 0,  v(5,0,t)  - 0 


u(5,ytt)  » Ug(5,t) 

In  a procedure  completely  analogous  with  steady  boundary  layers, 
equations  (14)  and  (15)  may  be  written  In  terms  of  Integral  relations  as 


follows: 


V - If  I?  ^ 


where : 


r6 

6*  - (1  - ~)d 

o^  6 

r« 

e - I — (1  - — ) 


^y|y-o 

The  solution  technique  employed  In  this  work  resembles  a procedure  used 

(12) 

by  Telpel  to  solve  the  classical  unsteady  momentum  Integral  equation. 

To  extend  the  range  of  pressure  gradients  for  which  this  procedure  Is 

applicable,  the  present  work  utilizes  a fifth  order  polynomial,  similar 
(13) 

to  that  of  Dryden  to  obtain  an  additional  degree  of  freedom  In 
choosing  the  polynomial  coefficients.  Thus,  the  general  form  of  the  u- 
component  of  boundary  layer  velocity  was  assumed  as 

u(C,y,t)  = a^  + b^r  + c^r^  + d^r^  + e^r^  + h^r^ 

+ Aa  cos  mt  [a^^  + b^^r  + C^^r^  + d^^r^  + e^^r^  + h^^r^] 

2 3 4 5 

+ Aa  sin  mt  ^ll’^  ^ll’^  ^ 

where  r ••  y/6(C,t) 

The  coefficients  In  equation  (21)  are  determined  by  matching 
the  following  conditions  on  the  wall  and  at  the  outer  edge  of  the  boun- 


dary layer: 


y - 0 : u - 0 and  3^.  - (_  + ^ — + qu^ 


'■7^‘ 
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The  scaling  function  q(t)  was  chosen  to  be  of  the  following  general  form: 
q(t)  ■ 1 + Aa  qj^(t)  - 1 + Ao  cos  ut  + C2  Aa  sin  ut. 

In  this  form,  and  C2  are  to  be  determined  such  that,  to  order  Aa, 
the  separation  point  always  occurs  at  a constant  value  of  Also,  an 
unsteady  shape  factor,  X,  Is  defined  to  be: 


, 3u.  _•  3u,  3u. 


The  unsteady  boundary  layer  thickness,  5(Ctt),  was  assumed  to  be  of  the 
form 


6(5»t)  • + 6j^Q(5)Aa  cos  ut  + 6j^j^(?)Aa  sin  ut 


where  5^(0  is  given  by  the  steady-state  analysis  and  and 


are  determined  from  the  unsteady  solution.  Substituting  gives  the  shape 
factor  In  expanded  form  as: 


X X + X-_  Aa  cos  ojt  + X-,  Aa  sin  u»t 
o XU  11 


du. 


du. 


d^u. 


du. 


10  o dC  o u,  o 1 dC'^  o 10  dC 

^0 

d^u 


“ ®o  dr""  - 

11  o d^  o u.  o 2^  d5^ 


- + 26  6 

o 11  dC 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


Recall  that  u,  , u.  , and  u,  were  given  as  functions  of  x only;  however, 
^o  ^10  ^11 

the  problem  solution  occurs  In  the  (^,y,t)  coordinate  system.  To  obtain 


values  of  any  component  of  the  potential  flow,  u^,  corresponding  to  a 
particular  the  Increment  form  of  Taylor's  Series  Is  applied.  For  exaiiq>le. 


consider  the  expansion  for  u^  (x)  given  below: 

o 


“6  " '*6  “ “6  ‘ 


= u^  (O  - Aa  5 qj^  0 

o d5 

Matching  the  boundary  conditions  gives  the  following  form  for  the  velocity 
profile. 

0 

where 

P(bj^,r)  - b^r  + (10  - 6bj)r^  + (8b^  - 15)r^  + (6-  3bj^)r^ 

G(b2,r)  - b2r  - j r^  + (|  - 6b2)r^  + (8b2  - |)r^  + (|  - 882) r^ 

In  equation  (32),  the  coefficient  of  r,  bj^  + 182*  may  be  arbitrarily 

chosen  so  as  to  Increase  the  range  of  X for  which  the  method  Is  valid. 

o 

Dryden^^^^  chose  bj^-  1.89  and  82  “ 0.11  for  his  analysis.  However,  for 
the  present  work  values  of  bj^  • 17/9  and  82  “ 1/7  were  chosen  because, 
when  compared  with  Dryden's  choices,  the  resulting  steady  lift  curve 
showed  less  deviation  from  that  obtained  by  the  classical  method,  yet 
Increasing  the  valid  range  of  X^.  Thus  u/u^  may  be  written  as 

— - F(r)  + X G(r)  + X,_G(r)  Ao  cos  ut  + X,,  G(r)  Aa  sin  ut 
Ug  o lU  11 


where 


17  4 3^14^15 

F(r)  - j- r - j r + j r 


o/^  1 12^9  3 5 4.1, 5 

G(r)  - yr-2r  +y^r  -y^r  +y^r 


f 


11 


The  model  for  the  boundary  layer  velocity  can  accommodate  time  varying 
phase  shifts  and  reverse  flow. 

Using  the  velocity  profile  (35),  the  Integral  relations  may  be 
expanded  as  follows: 


6*  ■ 6*  + 6*  Ao  cos  cot  + 6*  Aa  sin  cot 
o 10  11 


e - e + 6,-  Aa  cos  cot  + e,,  Ao  sin  cot 
o 10  11 


T ■ T + T Aa  cos  cot  + T Aa  sin  cot 
'^o  ''lO  ''ll 


(36) 

(37) 

(38) 


where 


0/17  . 0. 

'w 

o 0 


^o/lO  . /17  . 


®10  ^10, 


''lO  ^0  ' ^ 


1 “a  « 

. 0.11  . ,17  . V/_!ll  fll^i 


Substituting  Into  equation  (17)  and  discarding  terms  of  order  (Aa)  or 
larger  yields  the  following  set  of  governing  equations: 


order  (Aa)  : 


(200-^  V dT-iir 

U Oq 


(39) 


'^5 


168  dC 


3 

^0  5419  4 

^ X 1128X  8 ^ 

Q?  - 20^  + 6o  6o  (7  + (693)  (2520)  ''■  (693)  (49) 


fMl.  + 1.  X )—  —±—u  63^ 

231  ^520  49  o^d5  (693)  (49)  60  0 

3 

2 60  d0 

Q - — + — 2- 

168  u.  dC 

60 

n - 21233  _ , _ 4 _ ,2 

S u,  ^315  (1386) (1260)  0 (693) (49)  0 

0 


«0  d«0  443  . 4_  , . __4  ___  „ fi2 

■ ^ “60  “d?  49  (693)  (49)  "do  0 dC  ^ 


°0  0^443  4 ^ ^ 

^10  “ 693  ^2^  ■*■  T9 

Note  that  each  of  the  Q.’s  is  a function  of  the  steady-state  solution  only 
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In  the  governing  equations,  (40)  and  (41),  values  for  and  C2  are 

undetermined.  In  order  for  q(t)  to  properly  adjust  the  x-coordlnate  and 

C2  must  be  calculated  from  two  auxiliary  relations.  These  relations  are 

obtained  by  examining  the  unsteady  shear  given  by  equation  (38).  The  steady 

separation  Is  defined  as  the  point  for  which  t vanishes  or  equivalently 

o 

where  A * - 7(17)/9.  Thus  for  total  t to  vanish  at  this  point  for  all  time, 
0 w 

T and  T must  also  simultaneously  vanish  at  this  point.  This  occurs 

''10  ”11 

when  both  A^^q  and  Aj^j^  are  zero  at  the  separation  point.  Thus  equations  (29) 
and  (30)  give: 


du. 


10  ^ ,2,..  _!ii  _ 


'll(^s> 


du,  u,  d^u.  du, 

-d^  - ^ “ 60C25  ^ + 26o6n  ^ 

60 


From  equations  (42)  and  (43),  and  C2  may  be  calculated  such  that  q(t)  prop- 
erly scales  the  physical  x-coordlnate. 


CIRCULATION  CRITERION 

Determination  of  the  appropriate  potential  flow  about  the  oscillating 

body  Involves  specifying  correct  values  of  the  unsteady  circulation  components 

In  equation  (6) . The  criterion  used  In  the  present  study  to  determine  values 

(8) 

of  r^,  and  was  advanced  by  Sears''  • Sears  points  out  that  his 

(14) 

criterion  Is  an  extension  of  the  earlier  work  by  Howarth  based  on  a 
V..  ,,-..,-.(15) 


theorem  developed  by  Taylor 


The  Taylor-Howarth  Criterion  for  determining 
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the  proper  circulation  was  the  condition  that  the  net  vorticity  shed  into  the 
wake  from  both  the  upper  and  lower  surfaces  must  be  zero.  Sears'  modification 
of  this  requirement  to  account  for  unsteady  effects,  known  as  the  revised 
Taylor-Howarth  criterion,  is  given  by: 

dr 

* rate  of  vorticity  flux  into  the  wake 
Allowing  for  the  movement  of  the  separation  points,  the  criterion  in  (44) 


becomes 

dr  1-2  1 2 ^ ■'is  - ■‘*8 
« ■ 2 “S  2 Sd  + 28  ^ - “5  jT 

where  the  upper  and  lower  bars  denote  conditions  at  separation  on  the 
upper  and  lower  surfaces  of  the  body,  respectively.  Changing  the  criterion 
in  equation  (45)  to  the  (C,y,t)  coordinate  system  yields  the  following  three 
restrictions  on  Tq,  and  when  coefficients  of  1,  Ao  cos  wt,  and 

Aa  sin  ut  are  equated: 


“6o  " -6o 


„ „ - - . T ° “C-  C uu. 

T,  ■ (dF-,  - u-  u,  + C C,  u,  — 7p—  _ 7 2 8 6o 

1 11  6o  fiio  8 1 6o  dC 


+ u u - 5 C,  u 

-f.  -8  -1  - 

6o  filO 


<So 

6o  d? 


+ C,  £ (ou,  “ 0 
-2  ^8  -6o 


H H C,  u 

6o  6ii  6o 


-Si  58“5s„ 


From  equation  (46) , is  determined  by  Iteration  until  the  velocity  at  the 
upper  and  lower  separation  points  have  equal  magnitudes.  To  obtain  correct 


values  for  the  unsteady  components  of  circulation.  Iteration  Is  performed 
simultaneously  on  both  and  until  equations  (47)  and  (48)  are  satisfied. 

SOLUTION  PROCEDURE 

To  solve  for  the  steady  solution,  let  Z * 6^,  then  equation  (39)  may 

0 

be  written  In  the  following  form; 

/A  \ d^u. 

oq 

where  g(X  ) and  h(A  ) are  as  follows: 

0 o 


8(X^) 


h(X^) 


- 3.7777  + .81690  X - .1292  X^  - .00023559  x3 

O 0 0 

- .120097  + .00076101  X + .00029449  X^ 

0 0 

- .00050734  - .00023559  X 

0 

- .120097  + .00076101  X + .00029449  X^ 

0 0 


At  the  stagnation  point  u.  Is  zero,  thus  requiring  that  the  numerator  of  g(X  ) 

0 0 
0 

also  vanish  here.  The  physically  realistic  value  of  X^  for  which  g(X^)  becomes 

zero  Is  X^  - 5.0684.  Using  the  Initial  condition  obtained  from  applying 

1’  Hospital's  rule,  equation  (49)  is  Integrated  from  the  stagnation  point  to 

the  point  on  each  surface  where  t first  vanishes.  The  corresponding  values  of 

o 

u,  were  determined  and  the  criterion  of  equation  (46)  was  checked  to  decide  If 
0 

0 
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the  proper  value  of  had  been  chosen.  The  process  is  repeated  until  the 


proper  value  of  is  found. 

The  unsteady  solution  is  governed  by  the  linear  differential  equations, 

(40)  and  (41),  which  contain  the  unknown  constants  Cj^  and  C2  in  addition  to 

u,  and  Uf 

610  fill 

Correct  values  of  Fin  and  Fii,  which  fully  specify  u,  and  u,  along 

filO  fill 

the  body  surface,  must  be  determined  iteratively.  The  linearity  of  the  governing 
equations  permits  each  of  the  unknowns,  5]o  and  to  be  expressed  as  a sum 

of  three  coiiq>onents  as  follows: 

610  = ^fiio  + ^fiio  + C2  ^fiio 

fill  “ ^fill  + ^fill  + C2  ^611 

Substitution  for  610  and  fin  into  equations  (40)  and  (41)  permits  the  pair  of 
governing  equations  to  be  split  into  three  sets  of  two  equations  each  by 
equating  coefficients  of  1,  and  C2.  These  new  governing  equations  become: 


^filO  , , ““fill 

-5^  + ^fiio  Q2  + ^fiii  Q3  - — dT 


dC  ^6  de 


+ 

“filO 

“d^ 

1 

JO 

0 

d^fill 

du, 

filO 

d5 

fill 

*<2  ° 1 0 ^3  “ 

“ de 

^6 

(50a) 


dC  ^7 


(50b) 
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Qi  + ^<Sio  Qo  + Q,  + Qa 


(51a) 


-5^  Ql  + ^«ll  Q2  - ^«10  Q3  - Q5  = 0 (51b) 

dSfiio 

Ql  + ^«10  Q2  + ^^11  Q3  + Q5  = 0 (52a) 

-J^Ql  + ^611  Q2  - '610  Q3  + Q4  " 0 (52b) 

Note  that  equations  (51)  and  (52)  contain  only  Input  from  the  steady  solution 


(52b) 


and  may  be  solved  Independently  of  the  Iterative  process  necessary  to  determine 
Tio  and  Fu.  Thus  only  equations  (50),  which  are  solved  last,  depend  on  the 
choices  of  Fjo  and  Fji. 

Since  Qj^  vanishes  at  the  stagnation  point,  the  Initial  values  of  and 

^611  may  be  calculated  by  solving  the  two  algebraic  equations  obtained  by 
requiring  the  last  three  terms  of  equations  (51a)  and  (51b)  to  be  zero  at  the 
stagnation  point.  Similarly,  Initial  values  of  ^610  ^<511  are  obtained 

from  equations  (52a)  and  (52b) . Initial  values  for  the  derivatives  of  these 
four  unknowns  are  obtained  by  applying  1 'Hospital’s  rule  to  each  pair  of 
equations  and  solving  the  resulting  algebraic  relations.  Equations  (51)  and 
(52)  are  Integrated  by  Runge-Kutta  techniques  from  the  stagnation  point  around 
the  upper  and  lower  surfaces  until  the  points  of  steady-state  separation  are 
reached.  Values  of  the  variables  at  the  separation  points  are  determined 
and  retained  for  later  use  when  calculating  the  scaling  coefficients,  and  C2, 
on  each  surface. 


The  Iterative  portion  of  the  unsteady  solution  must  be  solved  for  each  set 
of  Fjo  and  F^j.  Once  again  Initial  values  for  ^610  and  ^611  are  needed  In  order 
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to  start  the  Integration  process.  However,  the  above  procedure  of  simply 

applying  1' Hospital's  rule  and  then  solving  the  resulting  simultaneous  equations 

Is  not  possible.  The  physical  reality  of  a moving  stagnation  point  (l.e.,  u. 

*10 

and  u.  may  be  non-zero  at  the  point  on  the  body  for  which  u.  vanishes)  prevents 
* 1 1 *0 

use  of  equations  (SO)  to  calculate  Initial  values  at  the  stagnation  point.  Close 

examination  of  equations  (50)  reveals  that  both  Q„  and  are  unbounded  when  u. 

o y on 


Is  zero.  Consequently,  unless  the  prescribed  potential  flow  components,  u. 


‘10 


and  also  vanished  at  the  steady-state  stagnation  point,  I'Hospltal's 

rule  would  not  be  applicable. 


In  the  present  study,  the  unsteady  stagnation  point  was  treated  In  a 
manner  similar  to  that  used  by  Rott^^^\  Once  a solution  was  obtained  for  the 
stagnation  region,  Initial  values  of  and  ^6^^  at  a small  distance  on 

either  side  of  the  stagnation  point  were  calculated  by  matching  the  shape  of 
velocity  profiles  resulting  from  Rott’s  procedure  with  those  used  in  the  un- 
steady momentum  Integral  technique.  Experience  showed  that  the  choice  of 
initial  values  for  ^6io  and  was  not  cirtical  because  boundary  layers  initially 

behave  In  an  almost  similar  fashion.  Thus,  any  reasonable  choice  of  starting 
values  for  ^<Sio  and  gives  Identical  results  after  Integrating  downstream 

to  the  separation  region.  Equations  (50)  are  Integrated  over  each  surface 
of  the  body  to  the  steady  separation  point  and  values  of  each 

separation  point  are  obtained  for  use  In  calculating  C^  and  C^  on  each  surface. 

To  solve  for  C^  and  C2,  £10  end  5xi  are  substituted  into  equations 
(42)  and  (43)  and  grouped  to  give: 


2 ^^10 

d^u. 

*0 

2 ^^10  ‘*”*0 

60  "W 

Cl  + 

60 

**“610  “611 


^610 


u 


60 


de 


(53) 


, 
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2 

du. 

do 

r*  M 

3a 

2 ®0 

d^u, 

do 

do 

dC 

«“^s 

do 

■ ^ dC^ 

^ *^“6 11  “fiio  2 


0 


(54) 


Once  values  of  Cj^  and  C2  are  obtained,  the  unsteady  circulation  criterion  Is 
used  to  decide  If  correct  values  were  chosen  for  Tio  and  Tij.  If  either  or  T2, 
equations  (47)  and  (48),  are  significantly  different  from  zero,  then  new  values 
of  Tio  and  Tn  are  chosen  and  the  solution  repeated  until  and  T2  are  both 
approximately  zero.  Prediction  of  new  values  for  Tjo  and  Tn  Is  rather  straight- 
forward after  the  results  of  and  T2  from  two  trial  values  are  obtained.  Since 

Tj^  Is  very  sensitive  to  changes  In  Tjo  and  insensitive  to  changes  In  Tji  and  T2 
exhibits  Just  the  opposite  behavior,  linear  extrapolation  Is  quite  adequate  for 
determining  new  trial  values  of  120  and  Ti}. 

Upon  completion  of  the  above  procedure,  the  unsteady  circulation  has 
been  determined  for  the  body  oscillating  at  a particular  ftequency  about  some 
mean  angle-of-attack. 

DISCUSSION  OF  RESULTS 

Once  the  components  of  total  circulation  have  been  determined,  the  two- 
dimensional  lift  coefficient  Is  then  obtained  from  the  following  expression 
relating  circulation  and  lift: 

Cj^(t)  - 2vr(t) 

A plot  of  steady  lift  coefficient  as  a function  of  angle-of-attack  was  formed 
by  repeatedly  solving  the  governing  steady-state  equations  for  various  and 
Is  presented  In  Figure  2.  A summary  of  steady-state  results  for  several  angles 


Table  1.  Summary  of  results  from  the  steady-state  solutions 


CIO 

To 

CD 

00 

0 

2 

0.0250 

0.8941 

-0.9502 

0 

6 

0.0695 

0.7399 

-0.9747 

0 

8 

0.0800 

0.5259 

-0.9836 

0 

9 

0.0750 

0.3427 

-0.9892 

of-attack  Is  presented  In  Table  1.  Corresponding  to  a particular  Table  1 
gives  the  steady  component  of  circulation  and  the  location  of  the  separation 
point  on  both  the  top  and  bottom  surfaces. 

Once  the  steady-state  results  are  available  as  Input,  the  unsteady 
governing  equations  are  solved  for  a nominal  angle-of-attack,  a , and  a | 

C t 

particular  frequency,  u.  Even  though  the  unsteady  governing  equations  are  j 

Independent  of  La,  Its  magnitude  Is  restricted  by  the  assumption  that  terms  of  | 

i 

second  order  or  higher  are  negligible.  Thus,  for  the  time  dependent  results  | 

presented  herein.  La  is  chosen  to  be  one-half  degree.  For  the  larger  values  J 

of  a^,  an  Increased  Aa  would  be  permissible. 

The  unsteady  results  for  a range  of  mean  angles-of-attack  and  frequencies 
are  presented  In  Table  2.  Notice  that  Table  2 Includes  the  scaling  coefficients, 
and  C2  In  equation  (23),  for  both  the  upper  and  lower  surfaces.  Also  the 
resulting  nonsteady  components  of  the  total  circulation,  r^o  Fu,  are  given 
for  each  and  frequency.  Using  the  value  of  F^  from  Table  1 and  values  of 
Fio  and  F^  from  Table  2 corresponding  to  a particular  a^  and  u,  the  unsteady 
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Table  2.  Summary  of  results  from  the  unsteady  solutions 


“o 

u 

"l 

£l 

^2 

^10 

^11 

0 

2 

0.10 

0.0666 

1.4501 

- .07900 

- .6665 

- .0231 

0.7209 

0 

6 

0.10 

- .3989 

5.4978 

- .07361 

- .4304 

- .0712 

0.4958 

0 

8 

0.01 

- .3372 

18.975 

- .00545 

- .4483 

- .0093 

- 0.0837 

0.10 

- 3.293 

18.446 

- .05350 

- .4541 

- .0899 

- 0.1018 

0.30 

- 8.264 

14.907 

- .14000 

- .4916 

- .2054 

- .2197 

0 

9 

0.01 

- 1.010 

44.466 

- .00033 

- .5325 

0.0061 

- 1.0952 

0.10 

- 9.752 

42.634 

- .00324 

- .5334 

0.0614 

- 1.0993 

0.30 

-22.85 

31.481 

- .00796 

- .5395 

0.1921 

- 1.126 

circulation  anf  lift  coefficients  may  be  easily  constructed.  In  Figure  3 lift 

hysteresis  loops  for  the  oscillating  elliptic  cylinder  corresponding  to 

several  nominal  angles-of-attack  are  presented.  The  effect  of  periodic  incidence 

changes  about  differing  values  of  is  apparent.  Note  that  for  * 2®,  the  loop 

is  very  thin  and  proceeds  in  a counterclockwise  direction.  The  results  at  6° 

are  very  similar  with  the  loop  progressing  in  the  same  direction  but  becoming  thicker 

as  the  region  near  , where  the  viscous  effects  are  more  prevalent,  is  approached. 

max 

For  a ■ 8°,  approximately  the  angle-of-attack  for  steady  C,  , the  loop  is  thick- 

max 

enlng  but  continues  to  exhibit  the  counterclockwise  direction.  In  the  steady- 
state  stall  region  at  ct^  ■■  9®,  the  significant  difference  is  that  the  loop  has 
changed  directions  and  now  proceeds  in  a clockwise  manner  making  the  direction  of 
the  hysteresis  loop  change  as  is  Increased  from  the  low  angles-of-attack.  This 
behavior,  however,  is  not  conmletely  unexpected  as  evidenced  by  the  experimental 
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(2) 

results  of  Lllva  et  al  and  presented  here  as  Figure  4.  The  results  presented 
In  Figure  4 were  obtained  from  two-dimensional  tests  of  a cylinder  of  airfoil 
cross  section  periodically  oscillating  with  an  angle-of-attack  Increment  of 
about  five  degrees.  Even  though  the  theoretical  results  of  Figure  3 are 
obtained  for  an  elliptic  cylinder  using  an  analysis  to  order  Aa  only,  a 
comparison  with  the  experimental  results  of  Figure  4 Is  enlightening.  Notice 
that  the  direction  of  the  experimental  hysteresis  loops  in  Figure  4 also  change 

from  counterclockwise  to  clockwise  when  a is  Increased,  as  do  those  of  Figure  3. 

0 

The  major  difference  however  is  where  this  direction  change  occurs.  In  the 
experimental  results,  the  loop  direction  change  is  indicated  by  a figure  eight 
behavior  in  the  hysteresis  loop  at  an  prior  to  that  for  steady  maximum  lift. 
However,  in  the  present  work,  the  loop  direction  changes  for  an  greater  than 
that  for  maximum  steady  lift.  This  difference  may  result  from  several  sources. 

One  of  these  is  revealed  by  examination  of  the  simplified  model  for  r(t)  in 
equation  (4) , which  Is  seen  to  be  Incapable  of  predicting  a figure  eight  type 
hysteresis  loop. 

Notice  in  Figure  3 that  the  predicted  hysteresis  loops  are  closely  aligned 
with  the  steady  lift  curve.  This  alignment,  which  Is  primarily  determined  by 
the  value  of  Fu,  seems  adequate  in  each  of  tue  results.  Of  more  Interest  perhaps 
is  the  behavior  and  Influence  of  Tig-  This  parameter  Is  largely  responsible  for 
the  thickness  of  the  hysteresis  loops  and  it  alone  determines  the  direction  in 
which  they  progress.  Positive  values  of  Fjo  produce  clockwise  directions,  whereas 
negative  values  of  Fig  produce  counterclockwise  loops.  In  Figure  5 the  effect 
of  frequency,  u,  on  Fjg  is  Illustrated  for  three  different  problems.  The  present 
work  is  labeled  "oscillating  cylinder."  The  curve  denoted  by  "oscillating  stream*^ 
was  obtained  by  applying  the  present  analysis  to  the  system  In  which  the  stream 


--  - 
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Is  perldocally  oscillating  about  a stationary  elliptic  cylinder.  Just  such 

a potential  flow  Is  obtained  from  equation  (2)  by  omitting  the  n term.  The 

z 

last  curve  was  obtained  by  casting  the  work  of  Moore^^^  In  the  notation. 

In  his  work,  Moore  also  analyzed  the  case  with  the  stream  oscillating  about  the 
stationary  elliptic  cylinder.  Judging  from  experimental  airfoil  results  In 
which  the  hysteresis  loops  collapse  at  the  large  frequencies,  If  test  data  were 
available  for  an  elliptic  cylinder,  the  general  slope  of  the  with  u curve 
would  likely  be  the  same  as  that  shown,  only  shifted  up  Into  the  quadrant  where 
Tio  is  positive.  Figure  5 also  Indicates  that  elimination  of  several  simpli- 
fying assumptions  made  by  Moore^^^  In  his  earlier  work  slgnflcantly  affects 
the  unsteady  lift  prediction. 

In  Figure  6 the  movement  of  the  separation  points  on  both  the  upper 
and  lower  surfaces  Is  depicted  for  * 9°  as  the  angle-of-attack  goes  through 
one  complete  cycle.  This  movement  Is  obtained  from  equations  (13)  and  (24) 
using  the  values  of  and  given  in  Table  2.  In  Figure  6 the  effect  of  In- 
creasing the  frequency  Is  to  produce  a type  of  phase  shift,  best  observed  at 
ujt  ■ ir,  in  the  movement  of  the  upper  separation  point.  Notice  that  the  lower 
separation  points  are  practically  stationary  whereas  in  Figure  6 the  upper 
separation  point  traverses  a segment  of  the  body  surface  equal  to  about  30 
percent  of  the  total  chord.  Also,  the  amplitude  of  the  upper  surface  separation 
point  movement  Is  vastly  different  for  each  half  cycle  of  the  angle-of-attack 
variation.  For  example,  with  Instantaneous  a(t)  ■ 9.5®,  the  forward  movement 
of  the  separation  point  Is  approximately  one  half  as  great  as  the  rearward 
movement  observed  when  a(t)  ■ 8.5®.  For  the  oi  ■ 0.10  curve  In  Figure  6 a 
type  of  phase  shift  Indicates  that  the  separation  point  lies  forward  of 
that  expected  from  steady  results.  Now  applying  the  logic  derived  from 
steady  boundary  layer  experience,  the  lift  is  expected  to  be  less  than 
the  corresponding  steady  lift.  That  Is  the  case  as  seen  by  examining  Figure  3 
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for  ••  9^.  Thus  the  movement  of  the  upper  separation  point  In  Figure  6 
appears  to  be  consistent  with  the  predicted  values  of  unsteady  lift. 

Perhaps  the  most  Interesting  aspect  of  this  analysis  Is  a comparison 

presented  In  Figure  7 of  the  present  results  with  the  case  where  the  stream 

Is  oscillating  about  a stationary  elliptic  cylinder.  For  values  of  of 

approximately  eight  degrees  or  less,  the  predicted  hysteresis  loops  for 

these  two  cases  are  very  similar.  However,  at  angles-of-attack  in  the  stall 

region.  Fig  Is  positive  for  the  oscillating  elliptic  cylinder  (clockwise 

direction  to  the  loop)  In  contrast  to  a negative  value  of  Fjo  (counterclockwise 

loop)  obtained  when  using  the  potential  flow  distribution  for  the  stationary 

elliptic  cylinder  In  an  oscillating  stream.  This  result  seems  to  support  the 

following  hypothesis  made  by  Moore^^^:  "It  may  be  that  different  directions 

of  hysteresis  should  be  expected  when  the  airfoil  oscillates  and  when,  as  In 

the  present  study,  the  stream  direction  oscillates."  Indeed  It  appears  that 

the  clockwise  direction  of  the  hysteresis  loop,  which  agrees  with  experimental 

airfoil  results.  Is  a direct  consequence  of  Including  the  term  in  the 

expression  for  u. (x,t)  which  alters  the  form  of  u,  to  produce  a positive  value 
0 Ojo 

of  Fig.  Physically  this  additional  term  seems  to  describe  the  Influence  on 
the  potential  flow  of  the  added  mass  effect  associated  with  an  accelerating 
body  In  unsteady  motion.  Thus,  the  experimentally  observed  phenomenon  that 
the  unsteady  lift  leads  the  periodic  angle-of-attack  change  may  be  at  least 
partially  explained  by  Including  the  term  In  the  potential  flow  description. 
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SUMMARY 

An  unsteady  momentum  Integral  technique  has  been  employed  to  analyze 
the  laminar,  two-dimensional,  unsteady  boundary  layer  on  an  elliptic  cylinder 
oscillating  periodically  In  a steady  freestream.  The  dynamic  lift  characteristics 
and  the  unsteady  behavior  of  the  boundary  layer  are  predicted  for  several  nominal 
angles-of-attack  over  a range  of  oscillation  frequencies. 

The  basic  approach  was  to  represent  the  boundary  layer  parameters  by  a 
steady-state  term  plus  an  unsteady  perturbation  proportional  to  Aa,  to  transform 
the  usual  boundary  layer  coordinate  system  such  that  the  location  of  the 
separation  points  appear  steady  In  the  new  coordinate  system,  and  to  relate 
the  rate  at  which  vortlclty  Is  shed  Into  the  wake  to  the  time  rate  of  change 
of  circulation. 

The  predicted  hysteresis  loop  direction  was  found  to  change  from 
counterclockwise  to  clockwise  as  the  mean  angle-of-attack  was  Increased  from 
two  degrees  to  within  the  stall  region.  This  direction  change  was  shown  to 
coiiq>are  favorably  with  experimental  results  for  cylinders  of  airfoil  cross 
section  with  the  principal  difference  being  that  the  experimental  results 
Indicate  the  direction  change  occurs  at  an  angle-of-attack  prior  to  that  for 
maximum  steady  lift,  whereas  the  present  study  predicts  a direction  change  at 
some  mean  angle-of-attack  after  maximum  steady  lift.  Probably  the  most  signifi- 
cant finding  resulted  from  a comparison  of  the  hysteresis  loop  predicted  at  a 
high  angle-of-attack  using  the  potential  flow  for  an  oscillating  body  in  a 
stationary  stream  to  the  loop  obtained  using  the  same  solution  techniques, 
only  changing  the  potential  flow  to  describe  a system  In  which  the  body  Is 
stationary  but  the  stream  oscillates.  In  each  case  the  hysteresis  loops  were 
similar  In  structure  but  with  opposite  directions.  The  clockwise  direction  of 
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the  predicted  loop  using  the  potential  flow  for  the  oscillating  body 
agrees  with  experimental  observations,  In  contrast  to  the  counterclockwise 
direction  obtained  for  the  system  with  the  body  stationary.  The  potential 
flows  for  the  two  systems  differ  only  by  a term  proportional  to  a(t)  which 
describes  the  Infliience  on  the  potential  flow  of  the  added  mass  effect 
normally  observed  when  a body  is  accelerating  In  unsteady  flow. 

An  obvious  shortcoming  In  the  present  model  Is  the  assumption  of 
laminar  flow,  particularly  on  the  upper  surface  of  the  elliptic  cylinder. 

By  considering  the  laminar  separation  point  as  the  point  of  wake  onset, 
an  excessively  large  portion  of  the  body  surface  Is  exposed  to  the  wake 
flow,  whose  effect  on  the  potential  flow  was  ^lssumed  negligible.  A logical 
extension  of  the  present  work  and  one  which  Is  being  undertaken  now.  Is 
to  Include  an  unsteady  turbulent  boundary  layer  model. 


28 


REFERENCES 


1.  Halfman,  R.  C.»  H.  C.  Johnson,  and  S.  M.  Haley.  "Evaluation  of  Hlgh- 
Angle-of-Attack  Aerodynamic-Derivative  Data  and  Stall-Flutter  Prediction 
Techniques."  NACA  TN  2533,  1951. 

2.  Lllva,  Jaan,  F.  J.  Davenport,  Lewis  Gray,  and  I.  C.  Walton.  "Two-Dimensional 
Tests  of  Airfoils  Oscillating  Near  Stall-Volume  1,  Summary  and  Evaluation 

of  Results."  USAAVLABS  Tech.  Report  68-13A,  U.  S.  Army  Aviation  Material 
Laboratories,  Fort  Eustls,  Virginia,  1968. 

3.  Ham,  N.  D.  "Aerodynamic  Loading  of  a Two-Dimensional  Airfoil  During 
Dynamic  Stall."  AIAA  J.  Vol.  6,  No.  10,  October  1968,  pp  1927-1934. 

4.  Ericsson,  L.  E.  and  J.  P.  Reding.  "D3mamlc  Stall  of  Helicopter  Blades." 
Paper  Presented  at  the  26th  Annual  National  Forum  of  the  American 
Helicopter  Society,  Washington,  D.  C. , 1970. 

5.  Crlml,  Peter.  "Dynamic  Stall."  AGARD  Report-AG-172,  1973. 

6.  Patay,  S.  A.,  "Leading  Edge  Separation  on  an  Airfoil  During  D3mamlc 
Stall."  ASRL  TR  156-1,  MIT  Aeroelastlc  and  Structures  Research  Laboratory, 
Cambridge,  Massachusetts,  1969. 

7.  Moore,  F.  K.  "Lift  Hysteresis  at  Stall  as  an  Unsteady  Boundary-Layer 
Phenomenon."  NACA  TN  3571, 

8.  Sears,  W.  R.  "Some  Recent  Developments  In  Airfoil  Theory."  J.  Aeronautical 
Sciences,  Vol.  23,  No.  5,  May  1956,  pp.  490-499. 

9.  Kochln,  N.  E. , I.  A.  Klbel,  and  N.  V.  Roze.  "Theoretical  Hydromechanics," 

John  Wiley  and  Sons,  Inc.,  New  York,  New  York,  1964,  pp.  320-332. 

10.  Moore,  F.  K.  "On  the  Separation  of  the  Unsteady  Laminar  Boundary  Layer." 
Boundary  Laver  Research,  Sprlnger-Verlog,  Berlin,  1958,  p.  296-310. 

11.  Rott,  N.  "Unsteady  Flow  In  the  Vicinity  of  a Stagnation  Point."  Quarterly 
Applied  Mathematics,  Vol.  13,  1956,  pp.  444-451. 

12.  Telpel,  Ingolf.  "Calculation  of  Unsteady  Laminar  Boundary  Layers  by  an 
Integral  Method."  Z.  Flugwlss,  Vol.  18,  1970,  pp.  58-65. 

13.  Dryden,  H.  L.  "Computation  of  the  Two-Dimensional  Flow  In  a Laminar  Boundary 
Layer."  NACA  TR  497,  1934. 

14.  Howarth,  L.  "The  Theoretical  Determination  of  the  Lift  Coefficient  for  a 
Thin  Elliptic  Cylinder.  Proceedings  of  Royal  Society  (London)  A,  Vol.  149, 
1935,  pp.  558-586. 

15.  Taylor,  G.  I.  "Note  on  the  Connection  Between  the  Lift  on  an  Airfoil 
In  a Wind  and  the  Circulation  Round  It.  Philosophical  Transactions  A, 

Vol.  225,  125,  pp.  238-245. 


29 


bi* 

Cl,  C2 


max 

L 

q(t) 

qi(t) 


Ti,  T2 


t 

u 

u 

u. 


«0 

“«10’  "fill 


w 

X' 

X 

X 

8 

Y' 


NOMENCLATURE 

major  axis  of  ellipse 
minor  axis  of  ellipse 

arbitrary  constants  chosen  to  model  the  velocity  profile 
scaling  constants  In  the  function  q(t) 
two-dimensional  lift  coefficient 
maximum  two-dimensional  lift  coefficient 
length  of  semi-chord  (ft) 

fiinctlon  of  time  used  to  scale  the  streamwlse  coordinate 

function  equal  Cj^  cos  u)t  + C2  sin  ut 

y/6 

test  fiuictlons  used  In  describing  the  unsteady  circulation 
criterion 

non-dimensional  time 
freestream  velocity  (ft/sec) 

non-dimensional  botindary  layer  velocity  in  the  x-direction 
potential  flow  velocity  (ft/sec) 
steady  component  of  the  potential  flow 
unsteady  components  of  the  potential  flow 
non-dimensional  boundary  layer  velocity  in  the  y-dlrection 
complex  velocity  potential 

coordinate  axis  coinciding  with  major  axis  of  ellipse 

non-dimensional  boundary  layer  coordinate  along  the  body 
surface 

x-coordlnate  at  separation 

coordinate  axis  coinciding  with  minor  axis  of  ellipse 

non-dimensional  boundary  layer  coordinate  normal  to  the 
body  surface 
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Subscripts: 


0 

10, 

11, 


angle-of-attack 
mean  angle-o£-attack 

fineness  ratio  of  the  elliptic  cylinder 
tmsteady  circulation  (ft^/sec) 
steady  component  of  the  circulation 
unsteady  components  of  the  circulation 
Increment  of  angle-of-attack 
boundary  layer  thickness 

form  the  boundary  layer  thickness  components  6io  and  dn 


displacement  thickness 

momentum  thickness 

boundary  layer  shape  factor 

transformed  streamwlse  coordinate 

location  of  separation  In  the  ^ coordinates 

wall  shear 

angular  velocity  (rad/sec) 
non-dlmenslonal  oscillation  frequency 
time  derivative 

overbars  denote  quantities  at  the  upper  separation  point 
underbars  denote  quantles  at  the  lower  separation  point 

steady  components 

unsteady  components  multiplying  La  cos  ut 
unsteady  components  multiplying  La  sin  ut 


-Mujjiuwim 
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a in  degrees 


teady  lift  hysteresis  for  several  values  of  wl 


Figure  4.  Unsteady  lift  hysteresis  for  a Vertol  23010-1.58  airfoil  at 
coos t^jot  oscillation  frequency  (Ref.  2) 
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variation  of  F,.  with  frequency  at 
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Figure  6.  Movement  of  the  upper  and  lower  surface  separation  points 

for  0-9° 
o 


